Charge Transport Through Open, Driven Two-Level Systems with Dissipation 
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We derive a Floquet-like formalism to calculate the stationary average current through an AC 
driven double quantum dot in presence of dissipation. The method allows us to take into account 
arbitrary coupling strengths both of a time-dependent field and a bosonic environment. We nu- 
merical evaluate a truncation scheme and compare with analytical, perturbative results such as the 
Tien-Gordon formula. 
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I. INTRODUCTION 

Coupled quantum systems with small effective Hilbert 
spaces are useful tools in order to study coherence, dis- 
sipation and the interaction properties of few-particle 
systems. In an electronic context, an example are cou- 
pled quantum dots 1 * 2 * 3 ^* 5 -, where strong interactions be- 
tween electrons 6,7,8 define a Coulomb blockade regime 
with tunnel-splitted many-body ground states separated 
from the remaining excited states. The ultimate limit 
of two states defines a two-level system for the charge 
degree of freedom, with electrons of a fixed spin tun- 
neling between two quantum dots. Studying transport 
and dissipation then leads to a non-equilibrium or 'open' 
(pseudo)spin-boson problem, where the coupling to ex- 
ternal reservoirs opens the path to investigate properties 
such as shot noised or decoherence in a controllable semi- 
conductor environment. 

Additional insight into the quantum dynamics of elec- 
trons can be gained by making the parameters of the 
problem time-dependent. When the time-dependence is 
slow, this can give rise to a variety of adiabatic phe- 
nomena such as charge p U mpingi2iiiii24i 3 Jii^4iSiiLiSii& ) 
adiabatic control of state vectors 20,21 , or operations rele- 
vant for quantum information processing in a condensed- 
matter setting 22, 23,24,25,26 i 27 i 28 . Different physics occurs 
in the high frequency regime where monochromatic time- 
variation induces photo-excitations, such as for coupling 
of AC fields to quantum doto 29 , 30 , 31 , 32 , 33 ! 34 , which has 
been tested experimentally 35 ! 36 ! 37 ! 38 ! 39 ! 40 recently. 

In general, AC driven systems^ 1 ^ and their applica- 
tion to various mesoscopic transporliMi 4 ^^ 4 ^^ 4 ^ 9 ^^ 
and tunneling 51 i 52 , 53 i 54 , 55 i 56 , 57 i 58 i 59 i 60 i 61 i 62 i 63 regimes 
have quite a long history, although the inclusion of 
interactions and correlations is a relatively new area. 
In low-dimensional systems, investigations have con- 
centrated on one-dimensional modela 6 ^ 5 ^* 6 ^^, the 
modification of Kondo-resonances by AC fields 6 ^ -, 
mean-field type approximations 7 ^ or exact studies of 
driven few-electron systems 7 ^ 3 .. 

In this paper, we combine AC driving with the dissipa- 
tive dynamics of a two-level system (double quantum dot) 



under transport conditions, i.e. in a situation where elec- 
trons in the Coulomb blockade regime can tunnel from 
reservoirs into and off two tunnel- coupled quantum dots, 
with the possibility to absorb from or emit bosons into a 
heat bath while simultaneously interacting with a classi- 
cal time-periodic electrical field. At first sight, combining 
such a multitude of possible interactions within one and 
the same model might look unsuitable for a useful the- 
oretical discussion. However, as we will demonstrate in 
this paper, it is possible to calculate experimentally rel- 
evant observables such as the time-averaged stationary 
current, with the help of the (heat bath) boson spectral 
density J{ui) as single, main input of the theory only. In 
particular, we show how within the polaron transforma- 
tion approach and for a given J (to), one can calculate the 
current for arbitrarily strong coupling to bosonic modes 
and an AC field. 

The paper is organized as follows: in section II, we 
describe the model Hamiltonian and derive a Floquet-like 
formalism for the stationary density operator. In section 
III, we compare analytical results for limiting cases with 
numerical data, and conclude with a short discussion and 
an outlook in section IV. 



II. MASTER EQUATION FORMALISM FOR AC 
DRIVEN DOUBLE QUANTUM DOTS 



In the following, we shall develop the general frame- 
work leading to explicit expressions for the stationary 
current through dissipative, driven double quantum dots. 
Our approach is in part similar to the treatment of closed, 
dissipative two-level systems with AC driving as reviewed 
by Grifoni and Hanggi— . Here, we generalise this ap- 
proach to take into account tunneling between the dots 
and the leads. In the non-dissipative case, this problem 
was treated by Gurvitz and PrageijZ 5 ^^ 6 . for non-driven 
double dots, and for coherently AC driven double dots 
by Stoof and Nazarcv^. 
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A. Model Hamiltonian 

We assume that the driven two-level system is de- 
fined in a double quantum dot device^. In the regime 
of strong Coulomb blockade, these can be tuned into 
a regime where the internal dynamics is governed by a 
time-dependent (pseudo) spin-boson model (dissipative 
two- level system^ 7 -), Hssit)- The latter describes one 
additional 'transport' electron which tunnels between a 
left (L) and a right (R) dot with time- dependent energy 
difference s(t) and inter-dot coupling T c (t), and is cou- 
pled to a dissipative bosonic bath (Hb — J2Q UJ Q a Q a Q)^ 



n SB (t) = 



2 



Q 

T c (t)a x +H 



9Q 
2 



(a-Q + aq) 



(1) 



The effective Hilbert space of double dot (without any 
coupling to electron leads or bosons) then consists of 
two (many-body) states \L) = \Nl + 1,Nr) and \R) = 
\Nl, Nr+1) and is defined by a pseudospin a z = \L)(L\ — 
\R)(R\ ee fiL-fiR and a x = \L)(R\ + \R)(L\=p + pl 

The effects of the bosonic bath are fully described as 
usual by a spectral density 



Q 



(2) 



where loq are the frequencies of the bosons and the qq 
denote interaction constants. When showing particular 
results we will be using 



J {(J) = 2auje- u/u 



(3) 



corresponding to a generic Ohmic bath. More realistic 
forms can be easily incorporated into our formalism but 
for simplicity, in this work we restrict ourselves to the 
Ohmic case. 

The coupling to external free electron reservoirs 

Hres — 

, £ fc Q c\ c k a is described by the usual tunnel 
Hamiltonian 



TlT 



E 



H.c), 



(4) 



s a = \0)(a\ (a=L,R). Here, a third state |0) = \N L ,N R ) 
describes an 'empty' DQD. Its presence leads to strong 
modifications both in the mathematical description as 
well in the physics of this problem, as compared to the 
case of an isolated spin-boson Hamiltonian. Here, the 
reservoir-related parameters of H re s and Ht have been 
assumed to be time-independent which again is an ap- 
proximation which might not be always fulfilled in exper- 
iments. Again, we concentrate on the simplest possible 
case in this work and neglect the effect of, e.g., a time- 
dependence in the external electro-chemical potentials. 
The full model as described by 



now offers the possibility to study non-equilibrium prop- 
erties of a time-dependent, 'open' dissipative two-level 
system. Note that in spite of the third, 'empty' state |0) 
we continue to use the term 'two-level system' here and 
in the following: although the presence of |0) leads to 
strong modifications of, e.g., the equations of motion of 
the density operator, it turns out that the internal dy- 
namics of the system is still closely related to that of the 
dissipative spin boson problem. 

The time-dependent spin-boson problem is in general 
characterised by the fact that bothe(t) and T c (t) are time- 
dependent. One can then investigate interesting effects 
such as adiabatic charge pumping, dissipative Landau- 
Zener tunneling^, or for the closed system (no coupling 
to the leads) the control of quantum superpositions 79 . 
Although this general time-dependence offers the richest 
spectrum of possible physical phenomena, one is clearly 
strongly restricted by the fact that nearly no analyti- 
cal solutions are available. In this paper, our goal is 
to develop a systematic theory for the stationary state 
of a somewhat simpler situation, i.e., the case where 
T c {t) = T c is constant, with the time-dependence solely 
contained in the bias s(t). 



B. Equations of Motion 

In the following, we treat the coupling to the reser- 
voirs within the Born and Markov approximation with 
respect to H i 33 i 80 , such that higher order effects like co- 
tunneling or the Kondo effect are not considered. This 
Born-Markov approximation becomes exact in the lim- 
iting case of infinite source-drain voltage 7 ^. Specifically, 
one sets the Fermi distributions for the left (right) reser- 
voir fi — 1 (/r = 0) whence the chemical potentials 
of the leads no longer play any role. Furthermore, the 
tunnel rates which are given by 



2tt 



\V k a fS(e-e ka ), a = L/R, 



(6) 



H(t) = HsB(t) + Hres + Ht 



(5) 



are assumed to be independent of energy. We mention 
that the generalisation to intermediate voltage regimes 
(finite bias) for double dots is a difficult and non-trivial 
problem even in the undriven case, which is why we only 
discuss the infinite-bias limit in this paper. 

The derivation of the equations of motion for the 
dot observables is now very similar to the non-driven 
case^S. The time-dependence of the Hamiltonian en- 
ters via the replacement of the phase factors e le( -* _t ' 
in the free un-driven time-evolution of the dots, by 
e if t , dse(s) f or ^ven case. Introducing the vectors 
A ee [fiL, nR,p,p'), r = r^ei (ei,...,e4 are unit vec- 
tors) and a time-dependent matrix memory kernel M, 
the equations of motion (EOM) can be formally written 
as& [(..) = n../>(t)], 

(A(t)) = (A(0))+ f dt' {M{t,t'){A{t')) +T) . (7) 
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This formulation is a useful starting point for, e.g., the 
calculation of shot noise. Note that in contrast to the 
undriven case, the memory kernel M depends on both 



times t and t' because there is no time-translation invari- 
ance in presence of driving. Explicitely, the equations for 
the dot expectation values read 



at 
o_ 

at 



n R ) t 



(P)t 



-iT c {{p) t - {pt) t } +T L [1 - (n L ) t ~ (n R ) t ] 
iT c {(p) t -tf) t }-T R {n R ) t 

dt l e itidse{s) 
dt l e -iSt,ds S (s) 



p<p) t / + iT c (n L ) t ,J C(t - f ) - iT c (n B ) t ,C*(t - t') 
(*Y^)t> ~ iT c (n L )A C*(t~t')+iT c (n R ) t ,C(t-t r ) 



(8) 



Here, half the decay rate (tunnel rate T R /2) of the system 
appears in the off-diagonal terms p and p\ acting as a 
source of dephasing due to tunneling of an electron out 
of the dot in either direction. Furthermore, the boson 
correlation function for a harmonic bath with spectral 
density J{uj), Eq. and at equilibrium temperature 
k B T =1/(3 enters, 



C(t) 

Q(t) 



-Q(t) 



(9) 



J(w) 



(1 — cos ujt) coth [ 



i sin tut 



In deriving, the equations for the off-diagonal elements 
(p(^), we used the polaron transformation (POL) and 
factorised the bosonic correlation functions from the dot 
operators in the equations of motions for the reduced 
density operator of the (pseudo) spin-boson system. This 
means that Eq. (JSJ is perturbative (though to infinite 
order) in the interdot coupling T c . 

Alternatively, one can perform a perturbation the- 
ory in the electron-boson coupling <?q, (weak coupling 
'PER' approach). In a calculation for an undriven dou- 
ble quantum dot, both approaches have been compared 
recently for the stationary currents^! and the frequency- 
dependent current noised. For the spin-boson problem 
with T R / L = 0, it is well-known that POL is equiva- 
lent to a double-path integral 'non-interacting blip ap- 
proximation' (NIBA) that works well for zero bias e = 
but for e 7^ does not coincide with PER at small cou- 
plings and very low temperatures. PER works in the 
correct bonding and anti-bonding eigenstate basis of the 
hybridized system, whereas the energy scale e in POL is 
that of the two isolated dots (T c = 0). This difference 
reflects the general dilemma of two-level-boson Hamilto- 
nians: either one is in the correct basis of the hybridized 
two- level system and perturbative in <?q, or one starts 
from the 'shifted oscillator' polaron picture that becomes 
correct for T c = 0. In fact, the polaron (NIBA) ap- 
proach does not coincide with standard damping theory 82 
because it does not incorporate the square-root, non- 



perturbative in T c hybridization form of the level splitting 
A = ^Je 2 + 4T C 2 . However, for large |e| > T c , A — > \e\, 
and POL and PER turn out to agree very well for the 
undriven case 



C. Stationary Quantities 

In a quantum system that is continuously driven by 
an external, time-dependent source, stationary quanti- 
ties can be defined for expectation values approaching 
a fixed point or a quasi-stationary, periodic motion for 
large times t. In particular, we will be interested in quan- 
tities like the time-averaged electronic current. It is then 
useful to split the time-dependent part off e(t) as 



e(t) = e + s(t), 



(10) 



and to introduce the Laplace transform f(z) = 
/ °° dte~ zt f(t) of a function f(t). The time-evolution of 
the isolated spin-boson system for T c = is governed 
by the correlation function C(t) = C*(—t). The Laplace 
transform of these, 



C E (z) = / dte- zt e let C{t) 
Jo 



c;(z) 



dte- zt e~ tet C*{t). 



(11) 



defines free propagators for quasiparticles in the uncou- 
pled dots and in absence of coupling to electron reser- 
voirs. In absence of electron-boson coupling, this simply 
describes the free time evolution of a particle described 
by the diagonal Hamiltonian etr z , whereas for non-zero 
boson coupling these become 'dressed' polarons. In addi- 
tion, the decay via the right reservoir at rate T R leads to 
a finite quasiparticle life-time and consequently a renor- 
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malisation of the propagators as 

C s (z) 



D e (z) = 

D*(z) = 



l + T R C E (z)/2 



E £ (z) E 

e;(z) 



l+T R C £ (z)/2 

C-e(z) 

i + r H £*(*)/2' 

(12) 



These expressions appear in the calculation in Appendix 
A, where Eq. JHJ is solved for the coherences (p) and (p^) 
in order to obtain two closed equations for the occupan- 
cies (n L/R ), 



zh L (z) - (n L ) 
zn R (z) - (n R ) 
K(z,t) 
G{z,t) 



dte 



dte 



(n L ) t K(z,t) - (n R ) t G(z,t) 



{n L ) t K(z, t) - (n R } t G(z, t) - T R n R (z) 



n L (z) - h R (z) 



dt'e- zt ' [T c (t + t')T*(t)D E (t') + T*(t + t')T c (t)D*(t')] 



dt'e~ zt [T c (t + t')T*(t)E e (t') + T*(t + t')T c (t)E*(t')} , 



where here and in the following we omit the (..) in the 
Laplace transformed expectation values to simplify the 
notation, and we defined 

T c {t) = T c e +l ^ ds ^ s \ T*(t) EET c e-^o^^ s ).(l4) 

Up to here the transformations have been valid for an 
arbitrary time-dependence in s(t). ^From now on, we 
specify to the time-periodic form 



e{t) = e{t + 2ir/n), 



(15) 



where 27r/f2 is the period of the time-dependent field (we 
further specify to a sinusoidal time-dependence of e(t) 
below) . 

We expect the system to approach an asymptotic 
quasi-stationary state. Then, the time-evolution of all 
quantities f(t) can be decomposed into Fourier series 



/(t)-/"(t) = £\ 



—inQ,t 



In 3 



(16) 



with multiples of the angular frequency f2 of the exter- 
nal field. Following Grifoni and HanggiZi, we decompose 



J 



K( Z ,t) = £#„,(*)< 
m 

G(z,t) = ^G m (z)e 



■tilt 



(13) 



K(z,t) and G(z,t) into Fourier series, 



(17) 



\asy 



The corresponding Fourier expansions {nL)t 
J2 m "™e- mint and (n R ) a t sy = £ m /i m e"™ m of the 
asymptotic occupancies can then easily be Laplace trans- 
formed, 



*r(*)=£' 



— (18) 

m 771 

and inserted back into Eq. (|13|) . Comparing the complex 
poles at z — —iMfl in the two equations for ul{z) and 
n R (z) and assuming that K m (z) and G m (z) are regular 
there, one obtains an infinite system of linear equations 
for the Fourier coefficients v m and p m , 



- iMVtv M = ~2_j WnK M -n(-iMn) - pL n G M -n{-iMVi)] + T L [S M ,o - v M - Pm] 

n 

[T R - iMSl] p AI = WnK M -n(-iM{l) - fi n G M -n{-iMQ.)} . 



(19) 



Upon adding these two equations, one has 

-l 



Mm 
= r M = 



1 



M ^ 0, (20) 



and Eq. (|19|l can be transformed into a single matrix 
equation for the coefficients v n , 



A 



(21) 



5 



where 



4 = 



& rn — 



(r L - irc.fi - r n T L ) 5 m , n 
K m ^ n (-imVL) + r n G m - n (-imtt) 

r L r R t r L 



G m (-imCl). (22) 



where the bar denotes the temporal average of the asymp- 
totic quantities over one period r = 2w/Q. This simple 
result means that the stationary current is determined 
by the Fourier component /io only. Note, however, that 
/Jo is part of the solution of an infinite set of the linear 
equations Eq. Using Eq. (fTT)|l and Eq. (J20J), one can 

express /io in terms of /in^o for the alternative expression 



D. Charge Current 

In the Master equation approach, the expectation val- 
ues of the electron current through the double dot is ob- 
tained in a fairly easy manner. One has to consider the 
average charge flowing through one of the three intersec- 
tions, i.e., left lead/left dot, left dot/right dot, and right 
dot/right lead. This gives rise to the three corresponding 
electron currents Jl(£), -ffl(i)) and the interdot current 
lLn(t)- iFrom the equations of motion, Eq. ©, one 
recognises that the temporal change of the occupancies 
( n L/R)t is due to the sum of an 'interdot' current cx T c 
and a 'lead-tunneling' part. Specifically, the current from 
left to right through the left (right) tunnel barrier is 

Ibit) = -eT L (n } t = -eT L [l - (n L } t - (n R } t } 
I R (t) = -eT R (n R ) u (23) 



and the interdot current is 

IlrV) = -<eT c {<p) t -<pt) t } 
d 
at 



(24) 



e-g t (n L ) t +lL(t). 



In the stationary case, all the three currents arc the same: 
adding the two equations Eq. Ijl9|l for M — 0, we first 
obtain 

Tfl/io = T L (1 - Mo - ^o)- (25) 

Using furthermore the Fourier expansion of {ni,/ R 
we recognize from Eq. (|23I24|I that 



asy 



I = I L (t) = I LR {t) = I R {t) = -eT R( x , 



(26) 



p *o(0) - £ n#Q [K_ n (0)/r n + G_„(0)] (i n 



R 



K o (0)/r o + G o (0) 



This form is in particular useful for the discussion of the 
Tien-Gordon limit below. 

E. Sinusoidal Time-Dependence 

In the following, we specify the time-dependence of the 
bias e(t) to a monochromatic sinusoidal field 



e(t) =e + Asin(Oi), 



(28) 



where the constant part is denoted by e. We introduce 
the notation, 



<7 e (-iw) = C e+a) (z = 0+) = C e+U 



(29) 



and correspondingly for the propagators D, D* , E, and 
E* , Ea. (|12|l . Then, invoking the decomposition of the 
phase factor into Bessel functions, 



T c (t + t')T*(t') = T C V^ + '' = (30) 

n) Jn ' (n 



and the definitions of the Fourier components K m (z) and 
G m (z), cf. Eq. (jT7|) and Eq. , one obtains 



K m (-im'n) = r m r c 2 ^ 

n 

G m (-im'n) = r m T^Y^ 



Jn 



Jn I ^ ) J 7i - 



J n 



J- n 



A 



D* 



e+(m' — n)Q T <Jn 



J„ I — I ' 



fi 



— (rn'+n)n 



(31) 



III. ANALYTICAL RESULTS tained, and then turn to a comparison with numerical 

calculations. 

In the following, we first discuss the limits where an- 
alytical results for the stationary current I can be ob- 
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A. Time-Independent Case 

For A = 0, i.e. in absence of the time-dependent 
(driving) part in e(t), we recover previous results^ for 
stationary transport in dissipative double quantum dots. 
One then has K(z,i) = K(z) and G(z,t) = G{z) 
such that K n (z) = G n {z) = for n ^ 0. Us- 
ing K (0) = 2Re[T c 2 4/(l + T R C E /2)], together with 
G o (0) = 2ReT 2 [C , * £ /(l +T R C E /2)}, after some algebra 
we re-derive the previous result 80 for the stationary cur- 
rent, 



B E = Re 




2Re(C £ 



l + T R C E /2\ 2 + 2T?B, 

C-e % 



(32) 



TrCs/2) 



1 



Fl 



(note the absence of the factor 2 in the definition of 
the rates here^S). The result Eq. (|32l) . which can be 
compared^ to an alternative derivation using perturba- 
tion theory in the boson coupling a, generalizes the case 
of elastic tunneling through double quantum dots to in- 
elastic tunneling with coupling to an arbitrary bosonic 
heat bath. For a = 0, we re-derive the Stoof-Nazarov ex- 
pression for the stationary current without dissipation, 



Iry — 



a=0 



^ 2 + ry4 + T c 2(2 + r iJ /r i ) 



(33) 



B. 



Lowest order T r : Tien- Gordon Result 



In the time-dependent case, we are able to derive an- 
alytical results by considering the limit of small interdot 
coupling T c , or large frequencies Q. These two limits do 
not yield identical results because apart from T c and f2, 
there are four other energy scales (bias e, rates T^, Tl, 
boson cut-off u> c ) in the problem. 

Considering Eq. Ij33(l for the undriven, non-dissipative 
current, lowest order perturbation theory in T c is valid for 
T cy f2 + r_R/r L < T R , \e\. The additional energy scale f2 
due to AC-driving requires that this condition is gener- 
alised to 



±0,1,2,.., (34) 



which indicates that at the resonance points e — nfl such 
a perturbation theory must break down, as is corrobo- 
rated by our numerical results discussed below. 

Considering the expression for hm in Eq. 119|) . one 
recognises that fij^ = 0(T 2 ) because the Fourier compo- 
nents of the functions K and G are proportional to T 2 , cf. 
Eq. i|13fl . Owing to the full expression Eq. (|27|l . the sta- 

TG 

tionary current in lowest order of T c is I = I + 0(T*) 
with 



■TG 



-eK o (0). 



(35) 



For a sinusoidal e(t) 
sion Eq. (|3*T)l yields 



e + A sin(Sli) , the explicit expres- 



jTG _ 



Re 



2C E j rn n 



1 



■#C, 



e+nVl 



(36) 



Note that Eq. (J2HJ) is the Tien-Gordon formula. This 
can be easily demonstrated by expanding the non-driven 
stationary current, Eq. I|32|l . to lowest order in T c , namely 
I = Io + 0(Tg ), such that, for the driven case: 



-TG 



iA=0 



(37) 



To lowest order in T c , the stationary current therefore 
is given by the Tien-Gordon formula: the current in the 
driven system is expressed by a sum over current contri- 
butions from side-bands e + nil, weighted with squares 
of Bessel functions. Note that the perturbative result 

I = —eKo(0), Eq. (J2EJl, does not refer to any specific 
form of the periodic function e(t); it is valid for arbitrary 
periodic driving when the corresponding Fourier compo- 
nent Kq(0) is used. 



C. Non-Adiabatic Approximation 

This approximation assumes that the frequency f2 is 
the largest energy scale in the problem, 



n^T c ,e,r R ,r L . 



(38) 



Onther.h.s. of the integral equation Eq. 11311 for n^/R^z), 

one then replaces the integral kernels K (z, t) and G(z, t) 
by their averages over one period of the AC field, 



k(z,t) 



2^ 



dtK(z,t) = K (z) (39) 



and similarily for G(z,t). The Fourier coefficients K n (z) 

and G n (z) with n ^ then vanish and one obtains / w 
jfast w h ere 



/-fast 



-eT R K (0) 



Tr + Go(0) + K o (0) [1 + Tr/T l } 



(40) 



We observe that within lowest order of the static tunnel- 
ing T c , Eq. (|40(l coincides with the Tien-Gordon expres- 
sion, Eq. (|35|l . which one obtains by setting Gq(0) oc T 2 
and A'o(O) oc T 2 to zero in the denominator of Eq. i|40[l . 
In fact, for the undriven case A = one can proved that 
the expression for the stationary current sums up an in- 
finite number of terms oc T 2 , a fact that can be traced 
back to the integral equation structure of the underly- 
ing master equation. Here, Eq. (|40|l demonstrates that a 
similar summation effectively can be achieved in the AC 
driven case. 
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D. Higher Order Corrections to Tien-Gordon 

In order to systematically go beyond the Tien-Gordon 
approximation, Eq. I|35|) . one has to perform an expan- 
sion of the current in powers of T 2 . This can be achieved 
by truncating the infinite set of linear equations, Eq. (|22|l . 
in order to obtain approximations for the n — 0, ±1, ±2- 
th side-band values of v n and fi n . The simplest way to 
do this in practice is by a numerical solution of these 
equations as discussed below. 

Barata and WreszinskiS have considered higher order 
corrections to dynamical localisation in a closed, coher- 
ent two-level system, i.e., without coupling to external 
electron reservoirs or dissipation. They found that the 
next order in perturbation theory given a contribution 
different from zero was the third order one, giving a con- 
tribution to a renormalisation of the tunnel coupling T c : 



= 



E 



2T C 3 
J2m + l (n) J 2 



(41) 



> S2 / 



J-2( ni +n 2 + l) (n) 



(2m + l)(2ra 2 + 1) 



We now recall our expression 
*o(0) = £ 



2RcL> F 



-nil 



(42) 



(and Go(0) correspondingly with D e+n ^ replaced by 
E e + n n), cf. Eq. (|31[) . which enter the Tien-Gordon result, 
Eq. (|35|1 . and the re-summed non-adiabatic approxima- 
tion Eq. (|4THl . We use the renormalised T c , Eq. l|4i"|l . in 
order to define a renormalised function (0) , 



K, 



(3) 



2ReD s+ „n, (43) 



and G (0) correspondingly. This yields an expression 
for the current, renormalised up to third order in T c , 
according to 



j(3) = 



-eT R K^{0) 



T R + G ( i ) (Q)+K { Q 3 >(0) [1 + T R /T 



(3), 



(44) 



In the following, we discuss and compare our above 
results. 



IV. DISCUSSION 

A. Comparison of Two Numerical Schemes 

In order to numerically solve the integro-differential 
system Eq. (JHJ, it is convenient to write 



cxi> I i J dse(s) 



(45) 




0.0001 



le-05 



FIG. 1: Average current through double dot in Coulomb 
blockade regime with bias e + A sin fit. Coupling to left and 
right leads Fl — Tr = T. Dotted lines indicate Tien-Gordon 
result, Eq.dSTt 



with (pt = et — (A/il) cos fit, remembering our choice 
e(t) — e + AsinSli. We then introduce the real and 
imaginary part of (p), use e lx = coscc+i sin a;, and specify 
to the Ohmic dissipation case for C(i), 



C(t) = \C{t)\e 

\c(t)\ = [i + M) 2 r Q 



2a arctanw c t 



(46) 



r(l + l//3w c + it/ (3) 



r(i + i/pu c ) 



4 a 



We have solved Eqs. JSJ numerically as a function of 
time, with the result for large times used to obtain the 
stationary current as a function of e. For each value 
of e, the time-dependent equations have been solved up 
to a fixed final time tf with a subsequent time-average 
over the interval [tf — At, tf]. tf has to be chosen suffi- 
ciently large, in particular for larger values of a. Conse- 
quently, one then also has to increase the number of steps 
to achieve sufficient accuracy of the data. We have used 
these numerical results to check our method for the sta- 
tionary quantities as obtained from truncating Eq. 1221) 
at a finite photo-sideband number, and found good agree- 
ment between both methods. Whereas the direct integra- 
tion of the equations of motion is somewhat slower than 
the truncation method, it has the advantage that it does 
not require analytic forms of the Laplace transform for 
the bosonic correlation functions C s , Eqs. (|11|) and Eqs. 
(|29H . The latter are required for the matrix scheme Eq. 
(12211 . In Appendix [B] we derive explicit expressions for 
zero temperature (T = 0) and Ohmic dissipation. Note 
that in contrast to usual l P(Ey theory, we require both 
the real part Re[C e (0)] = ttP(e) (where P(e) is the prob- 
ability for inelastic tunneling with energy transfer 3^), 
and the imaginary part of C e . 

In the following, we show numerical results obtained 
with the truncation method. 
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B. Photo-Sidebands (Coherent Case) 

1. Comparison with Tien-Gordon approximation 

In Fig. (JJJ, we compare the exact numerical result 
for the average stationary current with the Tien-Gordon 
expression, Eci. (|37[) . in the coherent case a = 0. One 
clearly recognises the symmetric photo-side peaks which, 
according to Eq. l|37|l . appear at ±nhe. The Tien-Gordon 
approximation overestimates the current close to these 
resonances, where terms of higher order in T c become 
important due to the non-linearity (in T c ) of the exact 
bonding and antibonding energies ± v / e 2 +4T C 2 of the 
isolated two-level system. This again confirms that the 
Tien-Gordon result is perturbative in the tunneling T c . 

2. RWA and Bloch-Siegert shift 



0.1 



a o.oi 



Vx.. x .. x ..xy \ 


A=0.3a 

A=0.6£2 

RWA 


\ \ 












T c = 0.25 £1 ^* 




r = 0.2 £2 





0.5 1 1.5 2 2.5 

elCl 

FIG. 2: Comparison between RWA, Eq. I|47jl . and exact 
result for first current side-peak. 



Close to the first side-peak, Stoof and Nazarov have 
used a Rotation Wave Approximation (RWA) to obtain 
analytical predictions for the first current side-peak. In 
this approximation, one transforms into an interaction 
picture where the fast-rotating terms with angular fre- 
quency ±f2 are transformed away, and terms with higher 
rotation frequencies (such as ±2f2) are neglected. The 
resulting expression for the current is 33 

A 2 T R (a 2 -4) 
SN c(cT 2 R .+ bA 2 ) w 2 + (e - e r ) 2 ' 1 ' 

with the resonance point er = -\/il 2 — 4T 2 and param- 
eters a = Q/T c , b = r^/r L + 2, c = a 2 + b - 4, and 
the half- width w = {a/[2^/a 2 - 4])^T 2 R + {b/c)A 2 . We 
compare JgN with the exact result in Fig.©. 

For smaller driving amplitude A, the agreement is very 
good but becomes worse with increasing A. The position 
of the side-peak resonance point, which is independent of 
A in the Stoof-Nazarov approximation Eq. 1)47(1 . starts 
to shift towards slightly larger values of the bias e. In 
fact, for stronger AC driving the RWA is known to break 
down: in isolated two-level systems, the first corrections 
to the RWA lead to the well-known Bloch-Siegert shift 84 
of the central resonance towards larger energies, which is 
consistent with the exact result in Fig. (J2J. 

C. Dynamical Localization and its Lifting 

In a quantum system driven by a periodic electric 
field, a phenomenon termed Coherent Destruction of 
Tunneling (CDT) (also denoted Dynamical Localization 
(DL)) occurs under certain parameters of the external 
fiel d 52 i 85 . The periodicity of the external field allows 
to write the solutions of the Schrodinger equation as: 
ip(t) = exp[— iejt]d)j(t) where ej is called the quasi- 
energy, and (j>j(t) is a function with the same period as 
the driving field: the Floquet state. 



When two quasi-energies approach degeneracy, the 
time-scale for tunneling between the states diverges, pro- 
ducing the phenomenon of coherent destruction of tun- 
neling (CDT)S. The time scale for localization is the 
inverse of the energy separation of the quasienergies. 

In the case of an isolated two level system driven by 
a monochromatic, sinusoidal field e(t) = e + Asin(Ot), 
Ea. (l28|) . CDT can be physically understood from the 
renormalization of the coupling T c of the two levels, 

T c -> T c , cff = T c J„ (A) . (48) 

This expression is obtained from first-order perturbation 
theory in the tunneling T^. At the first zero of the 
Bessel function Jo, namely when A/M7 = 2.4048..., the 
effective tunnel splitting vanishes, leading to a complete 
localization of the particle in the initial state. 

In the following, we discuss how stronger tunnel am- 
plitudes T c , the coupling to the external leads, and dis- 
sipation modify this picture. 

1. Current Suppression 

In Fig. ©, we show results for the average current 
and a = (no dissipation) in the dynamical locali- 
sation (DL) regime. Here, we define this regime by 
A = zo^l, where zq = 2.4048... is the first zero of the 
Bessel function Jo. For this specific value of the AC 
driving A, to lowest order in T c the average current is 
strongly suppressed for |e| < f2 as compared with the 
un-driven case A = 0. For small T c , this suppression 
is well-described by the Tien-Gordon expression (not 
shown here): since at A = zqQ, the n = term in the 
sum Eq. (|37|l is absent, the current is dominated by the 
shifted (un-driven) current contributions at bias e + nQ, 
with |n| > 1, which however are very small due to the 
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0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 

e/£l elQ. 



FIG. 3: Average current for AC driving amplitude A = 
zqQ. (zo first zero of Bessel function Jo) and various tunnel 
couplings T c . Coupling to left and right leads Yl — Yr = Y. 



resonance shape of the un-driven current. 



3 




L- ' 1 1 1 3 1 — 1 1 

-0.03 0.03 -0.03 0.03 

e/£2 e/£2 



FIG. 4: Central peak of average current through AC driven 
double quantum dot. Parameters T c = 0.1, A = zqQ (all 
rates in units of Q) LEFT: coherent case a = for different 
tunnel rates Y = Yl = Yr, dots indicate third order results 
Eq. 14411 . squares indicate the Tien-Gordon result Eq. 1371 for 
the case F = 0.005. RIGHT: disappearance of central peak 
with increasing dissipation a. 



2. Central Current Peak and Third-Order Result 

Surprisingly, however, the coherent suppression of the 
current is lifted again very close to e — 0, where a small 
and sharp peak appears. This peak becomes broader 
with increasing tunnel coupling T c , but its height is sup- 
pressed for increasing reservoir coupling T, cf. Fig. @ 
right. This feature is analysed in Fig. Q), where we 
show results for the central current peak around s = 
in the DL regime for coherent (a = 0, left) and inco- 
herent (a > 0, right) tunneling. As one recognises, the 
Tien-Gordon description (which is perturbative in the 
tunnel coupling T c ) breaks down close to e — where 
higher order terms in T c become important. As a matter 
of fact, for e = the only relevant energy scale of the 
isolated two- level systems is T c itself. In contrast, the 
third order approximation Eq. (|44|l reproduces very well 
the additional peak at e — 0, which indicates the impor- 
tance of higher order terms in that regime. At e = 0, the 
charge between the two dots is strongly de-localised in 
the undriven case, and this tunneling-induced quantum 
coherence persists into the strongly driven regime where 
its signature is a 'lifting' of the DL close to e = 0. 

The width of the corresponding current peak is deter- 
mined by the tunneling rate T. An increase of incoherent 
electron tunneling from the left lead therefore washes out 
the coherent lifting of the DL. This argument in empha- 
sised in the right part of Fig. (0} which shows that the 
central peak in the DL regime vanishes for increasing dis- 
sipation strength a. 




FIG. 5: Average current through double dot in Coulomb 
blockade regime with bias e + A sin Qt for various Ohmic dis- 
sipation strengths a at zero temperature. Driving amplitude 
A = fl for lines without symbols, A = zo^ {zo first zero of 
Bessel function Jo) for lines with symbols. Tunnel coupling 
between dots T c — 0.1Q, bath cutoff uj c = 500f2, and lead 
tunnel rates Yl = Yr = O.Olfl. 



D. Dissipation and Average Current 

1. Dissipative Photo-Sidebands 

As mentioned above, for simplicity we restrict our- 
selves to an Ohmic dissipative bath at zero temperature 
(T = 0) in this paper, leaving the finite temperature case 
or the case of more complicated spectral functions J(uj) 
for future work. 

For A = 0, we reproduce the analytical result Eq. i|32|) 
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FIG. 6: Average current through driven double dot for various 
AC driving amplitudes A and fixed dissipation a = 0.05, 
tunnel coupling T c = 0.1SI. 



and the corresponding inelastic current part for e > due 
to spontaneous boson emission&Sii. In Fig. (JjjJ, we show 
the stationary current as a function of bias e for various 
Ohmic dissipation strengths a at zero temperature and 
finite AC driving amplitudes A. For A = f2, apart from 
the central resonant tunneling peak, side-bands at e = 
ntt appear which reproduce the asymmetry of the central 
peak around e = 0. This asymmetry is a clear signature 
of the coupling to the dissipative environment strongly 
modifying the current even at zero temperature. 

The specific form of the inelastic current depends on 
the boson spectral density J(w)2!l. Note that in gen- 
eral, there is no monotonic dependence on the dissipa- 
tion strength a since the boson correlation function C e 
appears both in the denominator and the numerator of 
the expression for the current Eq. lol'l) . 



2. Dissipation and Dynamical Localisation in the Current 

If the AC driving amplitude A is increased towards 
ZqQ (zo is the first zero of the Bessel function Jo), one 
expects to enter the regime of dynamical localization and 
a strong suppression of the central current peak. In the 
coherent case a = (see above), resonant tunneling is 
usually strongly inhibited due to coherent suppression of 
tunneling. 

For a > 0, however, we find that the current sup- 
pression strongly depends on the static bias e: we find 
suppression for e > and, in general, larger values of 
the current for e < as compared to the case of smaller 
AC amplitudes A. We explain this feature in the follow- 
ing: the dependence of the average current on the driving 
amplitude A for fixed a is clearly visible in Fig. JBJ). A 
small driving amplitude A < 0.2 nearly does not change 
the current at all. However, the originally strongly asym- 
metric current curve becomes flattened out when A is 



tuned to larger values up to the dynamical localization 
value A = zqQ. There, the AC field nearly completely 
destroys the strong asymmetry between the spontaneous 
emission (e > 0) and the absorption side (e < 0) of the 
current. The central n = photo-band is completely 
suppressed and the dominant contribution to the current 
comes from the n = ±1 bands. For e < 0, the current 
for il > \e\ is due to photo-excitation of the electron into 
the first upper photo-sidebands and subsequent sponta- 
neous emission of bosons of energy E\ = Q — |e| to the 
bath. In contrast, for > e > 0, photon emission blocks 
the current because at T = there is no absorption of 
bosons from the bath. The remaining photon absorp- 
tion channel then leads to boson emission at an energy 
E% = SI + £, which is larger as compared to the case for 
e < 0, namely E% > Ex, and therefore has a smaller 
probability P(E) cx E 2a - 1 e- E/uJ % cf. Eq.JHU, leading 
to a smaller current. A similar argument can be used to 
explain why the current increases as one reduces e, say 
from e/n = to e/fl = —0.5. In particular, the strongest 
effect of the dissipative bath occur near one-photon reso- 
nance conditions, i.e. when e/Cl ~ ±1, where the current 
is regulated by the spectral function of the bosonic bath 
at very low frequencies, either by absorption (e < 0) or 
emission (e > 0) of a photon. These processes appear 
in the current as nonanalytic cusps reflecting the power 
law behavior of P(E). This has to be compared with the 
Lorentzian shape of the photo-sidebands in the absence 
of dissipation (Fig. 1). If one tunes to even larger values 
of A > zqVL, the central n = photo-band re-appears 
and the original strong asymmetry of the current curve 
is restored. 



V. CONCLUSION AND OUTLOOK 

Our results suggest that the combination of AC fields 
and dissipation in double quantum dots leads to a rich va- 
riety of non-trivial effects. In particular, we have shown 
that a time-dependent monochromatic field drastically 
modifies the dissipative inelastic stationary current, in 
particular for stronger AC driving in the dynamical lo- 
calisation regime. Corrections to the Tien- Gordon for- 
mula appear at larger tunnel coupling between the dots 
and become extremely important near zero bias in the 
DL regime, also in the non dissipative case. 

The method presented in this work has the benefit of 
accounting for an arbitrary dissipative environment via 
the correlation function C e . In the generic case, explicit 
analytical forms for this function are difficult to obtain 
and it might be easier to integrate the original equations 
of motion directly. Alternatively, one can numerically 
evaluate C e and use it as an input into our Floquet-like 
formalism. We also mention that the entire approach is 
based on the decoupling of the bosonic degrees of free- 
dom in the polaron transformed master equation. One is 
therefore always restricted to the range of validity of the 
NIBA (non-interacting blip approximation) of the origi- 
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nal spin-boson problemist. Discussing larger tempera- 
tures T should thus lead to more reliable results as com- 
pared to the 'test-models' C e which were discussed here 
for T = 0. 

A future extension of our approach should therefore 
be the derivation of a systematic perturbation theory in 
the electron-boson coupling, starting from the bonding- 
antibonding basis of the double dots. In a calculation 
for an undriven double quantum dot, such an approach 
has been successfully used recently to extract dephasing 
and relaxation times from the frequency dependent noise 
spectrum-,. 

Even for the coherent case a = 0, our results have 
shown that there are non-trivial effects due to the com- 
bined quantum coherence inherent in the double dot, and 
the coherence induced by the external driving field. In 
particular, we found systematic corrections to standard 
approximation such as the Tien-Gordon formula or the 
rotating wave approximation. The constituing quantities 
K m and G m of our theory, cf. Eq. I|27I31|) . describe dis- 
sipative tunneling cx T% of one additional quasi-particle 
between the two dots under the influence of the AC field, 
which again indicates that our approach is essentially 
perturbative in T c , although to infinite order and exact 
for a — 0. We showed that partial re-summations beyond 



the Tien-Gordon result are justified in a non-adiabatic, 
high-frequency approximation, but for the general case 
one has to rely on a systematic evaluation of Eq. $27). 
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SRC GR/R44690, the UK Quantum Circuits Network, 
the Nuffield foundation (T.B.), by the MCYT of Spain 
through the "Ramon y Cajal" program (R.A.) and 
project MAT2002-02465 (G.P. and R.A.). We acknowl- 
edge as well The EU Human Potential Programme under 
contract HPRN-CT-2000-00144. 



APPENDIX A: DOT OCCUPANCIES IN 
LAPLACE SPACE 

Here, we derive Eq. (|13[l for the occupancies {til/ r ). 
We define 

q(t) ee (p) t e-^o dse "W, gt(t) ee (p^te+'K dsi ^ 

(Al) 

This is inserted into the equations of motion in the time 
domain, Eq. JSJ, which upon Laplace transformation be- 
come 



zh L (z) - (n L ) = -i \ dte- zt {T c (t)q(t) - T*{t)q\t)} + T L 



n L (z) - n R {z) 



zh R (z) - (n R ) = i \ dte~ zt {T c (t)q{t) - T*(t)q^(t)} - T R n R (z) 



q{z) = — -q{z)C £ {z) - i 
q\z) = -^(z)C:(z) +l 



dt 



T c *(f ) (n L ) t ,C £ (z) - (n R ) t ,C*_ £ (z) 



dt'e- zt 'T c (t') [(n L ) t ,C;(z) - (n R ) t ,C^ e (z) 
I 



(A2) 



where we used the convolution theorem in the equations we obtain upon solving for q^> (z) and Laplace back- 
for q(z) and q*{z) and the definitions Eq. (fTTJl . Using transforming, 
the definitions for the propagators D and E, Eq. I|12|) 



q(t) = -i [ dt'T*{t')[{n L ) v D e (t-t')-{n R ) v E e {t-t')] 



q\t) = i I dt%(t') [(n L ) t ,D*(t - t') - (n R ) t ,E*(t - if)] , (A3) 

I 



involving the propagators in the time-domain. Insertion 
into Eq. I|A2|I yields 
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zh L (z) - (n L ) = - die zt 

Jo Jo 
ft 



dt'(n L ) t , [T c (t)T;(t')D e (t - + T* c (t)T c (t')D* £ (t - t')\ 



(A4) 



die 



-zl 



dt'(n R ) t , [T c {t)T* c {t')E e {t - t') + T*(t)T c (t')E*(t - t')} + T L - - h L (z) - n R (z) 



r 



At this point, it is useful to use a relation for a generalized 
convolution of a function K(t,t') and /(/;')> 



die 



dt'K(t,t')f(t') 



dte- zt f{t) / dt'e~ zt K(t + t',t) (A5) 
o Jo 

which can be easily proven by substitutions. Note that 
the usual Laplace convolution theorem is recovered from 
Eq. JSU if K (t, t') = K{t - t') is only a function of the 
difference of its two arguments. Eq. (|A4J| and a similar 
equation for nn(z) then lead to Eq. ((TBI . 



The second term is an analytic function of e. 
Now write 



~i\e\ 2a -\ e<0 

e 2a-l e -7ri(l/2+2a-l) £ > q 

e 2a ~ x (sin27ra + i cos27ra) s > 0. 



(B6) 



Recall the reflection formula for the Gamma function, 



APPENDIX B: CALCULATION OF THE BOSON 
CORRELATION FUNCTION 

Explicit expressions for the bosonic correlation func- 
tions C e , Eqs. (|ll|l and Eqs. I|29(l . which can be ob- 
tained in the zero temperature (T = 0) case for Ohmic 
dissipation. In this case, 



We have 

C(z) 



J {to) — 2aco exp(— lu/lo c ) 
C(t) = {l+iu c t)- 2a , i 



die~ zt {\ + iu) c t)~ 



2a. 



(Bl) 



(B2) 



{i( J j c )- 2a z 2a - 1 e- lz ' u "T{\ - 2a, -iz/w c ), 



where we used Gradstein-Ryshik 3.382.4 and T denotes 
the incomplete Gamma function. We set u) c = 1 for a 
moment to simplify notations and obtain 



C(-is) 



- i (-e) 2Q " 1 e- £ r(l-2a,-e) (B3) 



Note that e must have a small positive imaginary part 
here (Rez > in the definition of the Laplace transfor- 
mation): the incomplete Gamma function T(l — 2a, z) 
has a branch point at z = 0. However, we can use the 
series expansion 



r(l-2a,x) = T(l-2a) 



~ (_ 1 )n a .l-2a+n 

^ n!(l -2a + n) 

71 — K ' 



I -2a ^ 0,-1,-2,. 



to obtain 



(B4) 



C(-ie) = -«(-e) 

OO 



2c*-l -e 



r(l - 2a) 



(B5) 



^ n!(l -2a + n)' 

n=0 x ' 



2a^ 1,2,3,. 



Tfl-z) 



This yields 



£ > : C(-ie) 



r(z) sin 7rz 



_2a-l. 



(B7) 



r(2a) 



(B8) 



+ i 



Lr(2a) 



e 2a ~ x e~ e co\,2-Ka 



+ 6 E ^in\(l 



n=0 



(1 - 2a + n) 



e < : C(-is) = ie 



|2a— 1 



r(2a) sin2?ra 

OO 

^ n!(l 



n=0 



!(1 - 2a + n) 



From this, we can read off the real and the imaginary 
part of C(—ie). The real part is 



Red (-is) = irP(e) 



T(2a) 



e %a - 1 e- e 6(e). (B9) 



The imaginary part is 



Im(7(— ie) 



y t 

^ nl(l -2a + n) 

71 — V / 



TTISI 2 ^ 1 [ -1, 



(BIO) 



e < 



r(2a)sin27ra \ cos27ra, e > 
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